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Abstract 

Massive on-shell operator matrix elements and self-energy diagrams with outer gluon lines 
are calculated analytically at O(a^), using Mellin-Barnes integrals and representations 
through generalized hypergeometric functions. This method allows for a direct evaluation 
without decomposing the integrals using the integration-by-parts method. 



1 Introduction 



In the asymptotic region ^ m^, the heavy flavor contributions to the deeply inelastic struc- 
ture functions can be obtained from the corresponding massive operator matrix elements and 
the light flavor Wilson coefficients [1]. The massless Wilson coefficients for deeply inelastic scat- 
tering are known up to 3-loop order [2,3]. The heavy flavor contributions were calculated to 
next-to- leading order in [4] semi-analytically. A fast numerical implementation was given in [5] . 
Complete analytic results were derived only for the limit ^ m? for the structure function 
F^^^x.Q'^) to 0{al) [1] and ^(x, Q^) to 0{al) [6]. In both cases, the 0(af) massive oper- 
ator matrix elements are required. The asymptotic contributions cover all logarithmic and the 
constant terms, while contributions of 0{{m? / Q'^Y), k > 1 are not contained. In the case of 
the structure function F2{x, Q^), these terms yield a very good description already in the region 
> 20 GeV^, while for Fl{x, Q"^) this approximation only holds at large scales > 1000 GeV^. 
Since the heavy flavor contributions to the structure functions amount to 20-40 % in the small 
X region, cf. [7], and the scaling violations of these terms differ from that of the light parton 
contributions, their knowledge is essential for precision measurements of the QCD scale Aqcd in 
singlet analyses. 

In this letter we address a new compact calculation of the genuine 2-loop scalar integrals 
contributing to the massive operator matrix elements with outer gluon lines, based on the Mellin- 
Barnes technique [8-10] and using representations through generalized hypergeometric functions 
[11]. This approach allows to thoroughly avoid the use of the integration-by-parts method [12], 
which keeps the contributing number of terms low and yields very compact results. Moreover, 
we work in Mellin space to use the appropriate symmetry of the problem leading to further 
compactification. The complete calculation of the asymptotic heavy flavor Wilson coefficients 
will be presented elsewhere [13]. In the following, we will outline the principal method and 
present then the results for the seven contributing two-loop integrals in terms of nested harmonic 
sums [14, 15]. Some of the special sums needed are listed in the appendix. 

2 The Method 

The massive 2-loop diagrams considered are shown in Figure 1. The diagrams contain either 




Figure 1: Genuine 2-loop diagrams contributing to the massive operator matrix elements. 

three or four massive lines. The (S)-symbol in Figure 1 denotes the operator insertion of the 
corresponding local quark-gluon operators, see Figure 2. 

The diagrams can be decomposed into a *-product as described in Figure 3. We follow the 
calculation of Ref. [10], now generalized from massless self-energy diagrams to massive operator 
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Figure 2: QCD Feynman rules for the composite local operator insertions. 



matrix elements. 




Figure 3: The massive two-loop two-point diagrams with three and four fermion masses. The graphs are 
generated by inserting the three-point function, written as a double Mellin-Barnes integral, into the massive 
two-point function, i.e. performing the *-operation defined in the context of the Lie-algebra of Feynman 
diagrams [16]. The effect of this insertion is only given by a modification of the exponents of the two-point 
function, cf. [10]. 



Also in this case the above decomposition of diagrams can be achieved by applying the Mellin- 
Barnes representation ^ 

Let us consider Diagram I as an example. The corresponding gluing-product is depicted in 
Figure 4. Applying the Feynman-parameterization to the 2-point function yields 



^In his original contribution, Barnes notes that the contour integral representations (1) and those for more 
complicated integrands date back to Pincherle [17] Mellin [9] and Riemann [18], cf also [19]. 
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(27?)-^ {xikf + xiw? + 0:2(^1 — pY + a;2m' 

Here, A denotes a light-like vector with = 0, and D = A — 2e. Vi is the integer power of the 
respective propagator and i/jj... = z/j + z/j + . . . . The calculation is performed in the MS scheme 
and we factor out = exp[£(ln(47r) — 7^)] for each loop. ^ One shifts /ci — > /ci + X2P and the 
numerator term is decomposed as 

{A.h + X2A.pf = J2 ( J {A.hYix^A.pf'' . (3) 
1=0 ^ ^ 

All integrals over ki with {A.kiY, I > 1 vanish. 




Figure 4: Diagram insertion for graph I. 

One obtains 

^ ^ ^ r(z/4)r(z/i4 + AT - 1) ^ ^ 

In the next step, the three-point function is inserted in a similar way as in [10], since a part 
of the propagators are massive but no operator insertion occurs. The result for diagram I for 
Ui = l,Vi, can be expressed by a double Mellin-Barnes integral as 



iA7r)^{Aj)f^ 
(27r)2 r(l - 2e) 

X / da dr ^ ' r(-n A\ ^r(-r)r(r + 1) 

J ^i—ioo J ■y2—ioo v a -\- IJ 

r(a + r + l)T{a + r + 1 + e)T{-a - r - 2£)r(-(7 - r + e) 
^ r(a + r + 2)r(-a-r + l + A^) ' ^ 

The other six 2-loop integrals obey a similar representation. For fixed values of A^, one may 
calculate the Mellin-Barnes integrals using the mathematica-package MB by M. Czakon [20], 
which yields numerical values. They were given in Ref. [21], Table 1. ^ 
For the solution of one of the Mellin-Barnes integrals, relations like 

1 ^ T{a + s)T{b + s)T{d ~a-b- s)V{e - c + s)V{-s) 

3i^2[a, 6, c, e, 1] , 



^AU integrals are normalized to contain no mass-scale or factors of 27r in the final result. 
^Note that in [21] the spherical factor 5^ has not been factored out. 
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cf. [11], are used. The second integral is performed by applying the residue theorem. One obtains 



X 



r(s)r(i-s) 

+T{-e)T{l + e: 



r(l-2e) ^^r(^ + 2 + iV) 

r(j + 1 - 2e)T{j + l + e)T{k + j + 1 + N) 

r(j + 1 - e)r{j + 2 + N) r{k + j + 2) 

r(j + 1 + 2c)r{j + r(A- + j + 1 + c + iv) ' 

' r{j + l)r{j + 2 + e + N) r{k + j + 2 + e) 



(7) 



One first performs the £-expansion to the desired order. The infinite sums in (7) can be expressed 
through Mellin-type integral representations, partly using differential operators in the remaining 
summation and outer parameters. In some cases the starting values {k,j} = need a separate 
treatment. The corresponding Mcllin-integrals finally result into weighted harmonic sums. The 
calculations were coded in MAPLE. For fixed values of A^, simpler procedures are obtained which 
lead to analytic expressions for the expansion coefficients in e, cf. Table 2, Ref. [21]. ^ 

Except for diagrams VI, VII one may calculate the integrals in the above way for general 
values of N. 

Explicit representations for all diagrams could be derived using generalized hypergeometric 
functions for all diagrams. As an example, let us consider diagram VI, 



f dki f dk'i 



dk2 iAki)^-^-{Ak 



\N-1 



D 



{Aki-Ak2) 



x(47r)' 



2\l+2£ 



{ki — m'^){{ki—pY — m^){k\ — m^){{k2—p)'^—m'^Y{Jv2 — kiY 



(Ap)^-^r(l + 2£)(47r) 



r(i)r(j + 1 



2e 



27r 



N 




(8) 



B{l + e,l+j) 



r{j + 2 + 2e)r{j + l-e) 



3F2[l + 2£, -£,j+l;l,j + 2 + £; 1] 



After performing the expansion in e, the remaining sums can be carried out by suitable integral 
representations. Some of the sums required are listed in the appendix. 



3 Results 

All diagrams obey representations in terms of weighted harmonic sums in the present case. 
It even turns out that only single harmonic sums contribute. Not all the diagrams I-VII are 
independent. Due to the massless outer lines of the diagrams, graphs IV and V can be expressed 
in terms of diagrams I^ and II. Here I^ corresponds to the case where all scalar propagators have 
power one, i.e. 

IV = [i+(-in^xi„ (9) 

V = [1 + (-1)^^] J- X II (10) 
^In case of diagram and II checks could also be performed by nestedsums [22] at fixed A''. 
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For diagram III, even a closed expression for general values of e can be derived, 

^3 = 



dh f dk2 {47i)\Aki-Ak2)^-' ^-2U+2. 



(m 



(27r)^ J {27i)D {kf-m^){{ki-py-m^){{ki-k2y-m^)ki{k2-py 

2, Tt i-(-i)^ r(A^ + i + £)r(A^+ 1) 



(Ap)^-^r(l + 2£)(47r; 



sin(-7re) N{N + 1) T{N + 2 + 2e)T{N + 1 - e) 



(Ap)^-ir(l + 2£)(47r) 



2e 



N{N + iy 



1 2 

~s ^ N + 1 



+ 0(e] 



We summarize the results for the independent diagrams I-III,VI,VII in Table 1, expanding to 
0{6^) in the MS-scheme. The expressions in terms of harmonic sums were derived exploiting 
their algebraic relations [23]. 



Table 1: The analytic results for graphs I to VII for general values of N, with all Ui — 1, lb: 
ui = 2. 



la 

lb 

II 

III 

VI 

VII 



Sf{N) + 3S2(N) 

2N{N + 1) 
Si{N) ^ 3S2(N)/2 - Sl{N)/2 
N{N + 1){N + 2) 

Si{N) 



[l-( 



-1) 

.l)N 



N{N + iy 



1\ , g S^jN) 
e) N{N + IY 
1 

£ 



(iV+l)2(iV + 2) 
Sl{N)~S2{N) 



+ 



2N{N + 1) 



(iV + 1) 



4 
iV 



-1)^-1 



^i(iV) 



N 



2S,{N) 
m{N + l) ' N{N + 1)_ 



+ 



+ 



'iV2(iV+l)2 



+ 



Sf{N) - S2{N) + 25_2(iV) 2(3iV + l)Si{N) 



N{N + l) 



+ 



In terms of (single) harmonic sums 



N 



k=l 



the final expressions for the genuine 2-loop diagrams turn out to be extraordinarily simple. Due 
to the fact that we were not applying the integration-by-parts method, also the intermediary 
results remained rather compact. The results may be translated into x-space by inverse Mellin 
transformation using the Tables in [14]. On the other hand, one may work within the Mellin 
space representation continuing the expressions to complex values of N as described in [24] and 
construct the respective observables in analytic form. The inverse Mellin transformation is then 
performed by a single numeric integral. 

The diagrams are applied to express the unpolarized massive 2-loop operator matrix elements 
[1, 13]. As shown in [6, 13], nested harmonic sums will contribute in the physical result but they 
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are either due to one-loop insertions into one-loop diagrams or the four-propagator contributions 
to some of the topologies discussed here, which are obtained due to cancellation of numerator 
and denominator terms. 



4 Appendix 



We list some special sums which are typical for classes of sums to be derived for the present 
calculation. 



Si{k + z + N) 



1=1 



E 

i=l 



B{N + l,i) 

BiN,t) 
^ (iV + i + l)2 

(N + i) 



E 



E 



N-l 

E 



k=0 



N-l 



E 

k=0 



k 



L + 1 



■1)' 



k + 2 
k J{N-L + kf 



1 2Si{N) S,,,{N) (2^1 (iV) 2C3 
— — h 2- 



N m 

(2 - S2{N - 1) 

{-I f [2S.2{N) + (2 



2 N 



N 



N{N + 1) 

(2 - S^jN) 
N 



[2S_,{N) + C,] + 



N-l 

N{N + 1)^ 



+ (-1) 



AT 



C3 + S_s{N)-2 



S-2{N) 
N 



2N + 1 ^ N^ + 6N^ + 6N + 2 

;bi{N) + 



Here the symbol ai and Euler's Beta-function B{a, b) are defined by 

lim SUI 

r(a)r(6) 



(71 = lim SUN) 



B{a,h) 



r(a + h) 



(13) 
(14) 
(15) 
(16) 



(17) 
(18) 



Ar2(7V + l)2 ' Ar3(7V + l)3 

B{N -L,L + 2)[Si{N + L)-Si{N -L-1)] . (19) 



(20) 
(21) 
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